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Electric Polarization of Heteropolar Nanotubes 
as a Geometric Phase 
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The three-fold symmetry of planar boron nitride, the III-V analog to graphene, prohibits an electric 
polarization in its ground state, but this symmetry is broken when the sheet is wrapped to form a 
BN nanotube. We show that this leads to an electric polarization along the nanotube axis which 
is controlled by the quantum mechanical boundary conditions on its electronic states around the 
tube circumference. Thus the macroscopic dipole moment has an intrinsically nonlocal quantum 
mechanical origin from the wrapped dimension. We formulate this novel phenomenon using the 
Berry's phase approach and discuss its experimental consequences. 
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Physical properties of materials at the nanoscale can 
differ dramatically from their bulk counterparts. This is 
especially evident in the electronic properties, since the 
quantum behavior of electrons on this scale is sensitive 
to the size, shape and symmetry of the sample. Recent 
discovery of carbon nanotubes |L) provides a striking ex- 
ample, where metallic or semiconducting tubes of iden- 
tical compositions have only slightly different radii . 
Layered BN provides a III-V analog to these materials; 
it can be formed in single and multiwall nanotubes that 
have the same Bravais lattice as their graphene coun- 
terparts, but with inequivalent atomic species on its two 
sublattices H||. 

Here we show that the broken sublattice symmetry 
produces a macroscopic electric polarization in BN nan- 
otubes, dependent on their topology. Remarkably, this 
ground state polarization is an intrinsically nonlocal 
quantum effect that cannot be described by a classical 
theory. The sign and size of the longitudinal polariza- 
tion of the heteropolar tube are controlled by the bound- 
ary conditions on its electronic wave functions along its 
wrapped compact dimensions. We analyze this novel phe- 
nomenon by developing a quantum theory of the nan- 
otube polarization in terms of a geometric phase ■ 

A natural description of the electronic properties of 
this system is developed from an expansion of the tight 
binding Hamiltonian for the it electrons at small wavevec- 
tors q around the K and K' points at the corners of the 
Brillouin zone of the 2D hexagonal lattice p0pl| ] . Intro- 
ducing index a — ±1 for these points then leads to the 
long wavelength Hamiltonians Jl(] [l2|] 



H a (q, S, A) = ahv F q cr x + S a a y + A a z 



(1) 



where ct m are the 2x2 Pauli matrices 



q is the wavevec- 
tor along the tube axis, vf is the Fermi velocity and A is 
a site diagonal potential that distinguishes the B and N 
sites. The crucial parameter S a = ah~VFq± is a contribu- 
tion to the electronic gap that comes from the quantized 
crystal momentum q± along its azimuthal direction. 



We can interpret Eqn. ([!]) as the Hamiltonian for an 
effective spin 1/2 particle interacting with a "magnetic 
field" defined by the parameters (q, 5, A). H a can be 
diagonalized by the rotation U a (q, S, A) = exp(iCl a -er/2), 
which aligns the vector 

(2) 



(units hvp — 1) in the z direction. The valence band 
eigenfunctions v a of Eqn. (|lj), r elevant here, have ener- 
gies E a (q) = —\Jq 2 + 5 a + A 2 . They are obtained by 
the rotation of the down spinor vo — (0, 1) according to 
v a (q,6,A) = exp(ir2 Q • a/2) ■ Vo, and correspond to the 
spin down state quantized along the local Q a axis. 

The electric polarization of an extended system is ill 
defined as an intensive quantity, but its changes are well 
defined For an electronic Hamiltonian containing 

a control parameter H(X), the difference in polarization 
Ap between initial and final states at Xi and Xf can be 
obtained by integrating the differential changes in p as 
one adiabatically varies the control parameter A 



Ap 



dX dx 



(3) 



We evaluate Ap on the path from the C nanotube, which 
is nonpolar by symmetry, to the heteropolar BN nan- 
otube which can have a macroscopic dipole moment. We 
take as a control parameter the antisymmetric site diago- 
nal potential A, where Aj = in the C nanotube. There- 
fore, integrating Eqn. (^) over the range < A < Ay, 
gives an expression for the electric polarization of the 
heteropolar BN nanotube with the ionic strength Af. 

Equation (^) can be evaluated by studying the vari- 
ation of — i (v a \d/dq\v a ) as a function of A. Summing 
over a and wave vectors q, we find 
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The valence states v a adiabatically follow the direction of 
the vector — Q a , so that dv a /dq = (— i /2)(d(fl ■ a) / dq)v a 
and dv a /dA = {-i/2){d(Q ■ a)/dA)v a . Therefore, the 
polarization is 



p - = -y.r 

e 4tt ^ .L 



dA 



dq 



(v a \a y \v a ) 



S 2 + A 2 



(5) 



Equation (^|) has a simple geometric interpretation. 
The effective "spin" representation for each valence state 
v a defines a unit vector field directed inward along the 
radial direction at each point (q, <5, A) in the parame- 
ter space for 7i a . The two surface integrals in (^) give 
the flux of a radial inverse square field in this parame- 
ter space that links through the two rectangular loops, 
shown in Fig. [l]a. These are located at 5 a= ±\ = ±hvp q±, 
oriented parallel to the q — A plane and extending from 
-7r < q < n and < A < A/. 





FIG. 1. (a) The electric polarization is the flux of a vec- 
tor field linking through two rectangular loops in the space of 
Hamiltonian control parameters, (b) The flux through each 
rectangle is proportional to the solid angle subtended by the 
loop, which is also swept out by the vector Cl when trans- 
ported around the perimeter of the loop. 

This flux is related by Stokes' theorem to the line inte- 
gral of an effective vector potential {v a \d q ^\v a ) around 
the loop perimeters, shown in Fig. ^b. Thus, the elec- 
tric polarization is proportional to the solid angle that 
is "swept out" by the effective spin as it is adiabatically 
transported around the loops. 

Introducing an energy cutoff at the bandwidth W = 
Trhvp, the double integral in Eqn. (||) gives 



dA 



dq 



ad 



(q 2 + S 2 + A 2 ) 3 / 2 



= 2a arctan 



WA 



sVs 2 



w 2 



(6) 



It is intuitively clear that by interchanging the B and N 
sublattices, i.e. negating A, the sign of p should change. 
Much more profound is the fact the dipole moment in 
Eqn. (JsJ) is also odd in S, which thus reflects its essen- 
tial quantum mechanical origin "rooted in the wrapped 
dimension" . 

The macroscopic polarization is obtained by summing 
the allowed values of 6 for each occupied subband in 
Eqn. (Q). When a BN sheet is wrapped into a tube 
the transverse crystal momentum is quantized q_L, n > and 
this restricts the allowed values of 6 in the sum. The 
transverse momenta q±. n depend on the subband in- 
dex, n and the wrapping vector Cmn — MT\ + NT2 
of the BN nanotube, where T\ and T2 are the primitive 
translation vectors for the hexagonal lattice p| ^P^l ■ A 
tube with wrapping indices [M, N) has a chiral index 
v = mod (M — N, 3) and its transverse momenta are 
quantized to the values q_i_ t n — (27r/|C M iv|)(n+5Sgni/), 
so the gap parameter is 5 an = ahvF q±,n- Thus sum- 
ming Eqn. (H) over subbands finally gives the net electric 
polarization 



1 



V 



2tt ^ 



a arctan 



IFA 



S an y/S 2 an + A 2 + W 2 



(7) 



Note that the summand is odd in $cm, but even under 
the interchange a «-> — a. Thus contributions to the elec- 
tric polarization from states near the K and K' points 
are additive. However, a nonzero electric dipole moment 
occurs only when the momentum distribution around ei- 
ther point is odd in S an . Thus the quantization of the 
transverse crystal momentum is the critical feature that 
controls the electric polarization of the BN nanotube. 

Eqn. (0) predicts that BN nanotubes with crystal mo- 
mentum distributions even in q± are unpolarized. This 
happens in the limit of the flat BN sheet, for which the 
macroscopic dipole moment is obtained from an integral 
over all two-dimensional wavevectors, and for tubes with 
chiral indices v = 0, for which the kinematically allowed 
momenta occur in ±S an pairs. A simple example of the 
latter occurs for the "armchair" M = N structures with 
the mirror symmetry (or invariance under M <-> N) [ [l4[ . 
Interestingly, Eqn. (^) predicts a zero polarization even 
for "nonarmchair" tubes with v = 0. 

However, for tubes with v 7^ the i5 <-> — <5 symmetry 
is broken by the fractional quantization of the transverse 
crystal momentum, and consequently these tubes have a 
net electric polarization with the sign determined by v. 
This is a remarkable result, as it implies that two nearby 
structures (M, TV) and (M+l, N) can have opposing elec- 
tric dipole moments, even though their structures in the 
tangent plane of the tube are nearly identical. 

A striking illustration of this effect is given in Fig. ||, 
in which we plot the electric polarization (the dipole mo- 
ment per unit length in one dimension), calculated as a 
function of the symmetry breaking potential A, for three 
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different "zigzag" tubes [Q with the nearby wrapping 
indices (M, N) = (17,0), (18,0) and (19,0). Skeleton 
lattice structures for these three tubes are shown in the 
insets. We use the following parameters representative 
of BN: hv F = 5.4 eV A, W = 10 eV and A = 2.5 eV, and 
plot the polarization in units of the elementary charge e. 




A(eV) 



(19,0) 

FIG. 2. The polarizations (dipole per unit length) for three 
heteropolar zigag tubes are plotted as functions of their ionic- 
ity parameters A. The curves are labelled by the tube wrap- 
ping indices (M, N), and the insets give their skeleton lattice 
structures. The (18, 0) tube has no net polarization while 
its neighbors (18 ± 1,0) have nonzero electric polarizations of 
opposite sign. 

Figure || shows that the electric polarization is zero 
for the (18, 0) tube, and it is inverted for the (17, 0) and 
(19,0) tubes, despite the fact that the three tubes have 
nearly identical local atomic structures in their surface 
planes. This reflects the interesting fact that the dipole 
moment arises as a nonlocal quantum effect that is con- 
trolled by the continuity of the electronic wavefunctions 
around the tube circumference and not any local prop- 
erty of the BN lattice. In this way, the circumferential 
boundary conditions for the electronic Bloch states can 
ultimately distinguish the three structures by their chiral 
indices v. 

The polarization (when nonvanishing) for BN tubes 
corresponds to a dipole moment of nearly 0.5 D per unit 
cell. Its sign, determined by the chiral index is, alternates 
with a three sublattice period as a function of the wrap- 
ping indices M and N. This effect is the heteropolar ana- 
log of three sublattice symmetry that distinguishes con- 



ducting and semiconducting behavior of graphene nan- 
otubes Thus semiconducting graphene nanotubes 

with v ^ may be classified into two families that are 
distinguished by the sign of their gap parameters. This 
controls the sign of the polarization induced by adiabatic 
introduction of the ionicity parameter A in Eqn. (□). 

The results derived here from the Hamiltonian (pi) give 
the contribution to the electric polarization from the 7r 
derived electronic states. In general, one can expect ad- 
ditional contributions to the electric dipole moment from 
lower lying a bands. However, these contributions are not 
expected to show a strong dependence on the chiral index 
and should be smaller than the effects derived from the 
more polarizable 7r manifold. Thus the periodic depen- 
dence of the electric polarization on the wrapping vector 
arises from the phase matching of the 7r-derived Bloch 
states around the tube circumference. It will be useful 
to quantify the relative tt- and a- derived contributions 
to the electric polarization within a complete microsopic 
theory 

Experimental observation of the macroscopic dipole 
moment of a heteropolar tube is complicated by the pos- 
sibility that extrinsic or surface charges accumulated at 
the tube ends could mask the predicted intrinsic effects. 
Nevertheless, the polarization p in Eqn. defines (mod- 
ulo e) the surface charge that must appear at the ends 
of the v 5^ heteropolar nanotube, independent of the 
atomic structure of the termination. Thus two-thirds of 
heteropolar single wall nanotubes should contain a uni- 
versal but nonintegral surface charge, determined by the 
bulk electronic structure. Similarly junctions between in- 
equivalent tube segments are predicted to localize inter- 
face charges 0, ±p or ±2p modulo e. 

Multiwall heteropolar nanotubes would favor combina- 
tions which compensate the macroscopic dipole moment. 
For a double wall tube of outer radius a and inner radius 
6, one can estimate the electrostatic contribution to the 
surface energy 



U, = 



2z 2 e 2 
ir\a — b\ 



K 



/ Aab \ 
\(a-b) 2 ) 



(8) 



where ze is the effective bound charge accumulated at 
the tube end and K is the complete elliptic integral of 
the first kind. Taking z w 0.2, a w 10.4 A and b w 7.0 A 
gives U s ~ 65 meV. Note that the structure of the double 
wall tubes grown under experimental conditions can be 
influenced by various kinetic factors as well. 

Our model can be experimentally tested by measur- 
ing the change in polarization induced by varying control 
parameters in the Hamiltonian. Of particular interest is 
the effect of an elastic strain e linearly coupled to the 
gap parameter 5 — > S + Ae. Thus the heteropolar tube 
is a molecular piezoelectric with the response constant 
z* = dp/de = Xdp/dS given by 



3 



AW(26 2 + A 2 + W 2 ) 



2tt (5 2 + A 2 )(S 2 + W 2 )VS 2 + A 2 + W 2 



j-g-j of the important role of the geometric phase in quantum 
mechanics. 



The gap parameter S can be linearly coupled to long 
wavelength compression, extension or torsion of the tube 
depending on the wrapping indices Jig] . These strains are 
the tube analogs to the orthorhombic and shear strains 
of an isolated graphene sheet. 




FIG. 3. The piezoelectric constants for the heteropolar 
zigzag (17, 0),(18, 0) and (19,0) tubes are plotted as a func- 
tion of the ionicity parameter A. 

In Fig. H we plot this piezoelectric constant as a func- 
tion of the ionicity parameter A for the heteropolar 
(17, 0), (18, 0) and (19, 0) nanotubes. Although the static 
dipole moments p/e for the (17, 0) and (19, 0) tubes have 
opposite signs, their piezoelectric constants z* /e have the 
same signs. This is because the couplings of the gap to 
strain are nearly identical for the two structures. The 
piezoelectric constants should be large for the BN tubes, 
since we expect that A ~ 1 eV, similar to carbon nan- 
otubes. Note that the unstrained (18,0) tube has no 
static dipole moment but it has an exceptionally large 
piezoelectric constant, that diverges proportional to 1/A 
as A — ► 0. This occurs because A = marks a critical 
point, where the gap vanishes for tubes with chiral index 
v = 0. Figures |||^ also demonstrate that BN realizes the 
"strongly ionic" limit of this problem, in which the static 
dipole and its piezoelectric coefficient are essentially sat- 
urated. 

The electric dipole moments of heteropolar tubes pro- 
vide crucial information for understanding their struc- 
tures, elementary excitations and phase behavior p6| , pT{ . 
They are related, for example, with an assortment of 
new photogalvanic effects [jllj . It is inevitable that other 
physical phenomena and applications will derive from the 
unique properties of nanotubes with controllable electric 
polarization. This macroscopic phenomenon also pro- 
vides a beautiful and experimentally relevant illustration 
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